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Models of the spontaneous emission of photons coupled to the electronic 
states of quantum dots are important for understanding quantum in- 
teractions in dielectric media as applied to proposed solid-state quan- 
tum computers, single photon emitters, and single photon detectors. The 
characteristic lifetime of photon emission is traditionally modeled in the 
Weisskopf-Wigner approximation. Here we model the fully quantized 
spontaneous emission, including near field effects, of a photon from the 
excited state of a quantum dot beyond the Weisskopf-Wigner approxima- 
tion. We propose the use of discretized central-difference approximations 
to describe single photon states via single photon operators in 3+1 di- 
mensions. We further show herein that one can shift from the traditional 
description of electrodynamics and quantum electrodynamics, in terms of 
electric and magnetic fields, to one in terms of a photonic wave function 
and its operators using the Dirac equation for the propagation of single 
photons. 

The field of quantum computation (QC) and quantum infor- 
mation technology (QIT) has recently experienced escalated 
activity in the search for sources of photonic states coupled to 
their quantum sources [1,2]. The observed increase in activ- 
ity is in part due to the suggestion that quantum information 
processing based on the electron spins of quantum dots (QDs), 
coupled through optical modes of a micro-cavity, [1] could im- 
prove on schemes based on the energy states of trapped ions 
[3] and nuclear spins in chemical solution [4, 5]. Some ad- 
vantages of a QC scheme based on the suggested semiconduc- 
tor quantum dot arrays may include greater scalability, longer 
spin decoherence times, longer coherence lengths, and fast in- 
teractions mediated by photons [1]. In this contribution we 
present a model for describing the coupling between a single 
photonic state to the spin states of a quantum source (such as 
a QD) through optical modes present in a micro-cavity. Ad- 
ditional applications of this model may include the design of 
devices aimed at single photon emission [6], single photon de- 
tection [7, 8], quantum teleportation [9, 10], quantum comput- 
ing within a quantum network [2], and quantum cryptography 
[11, 12, 13, 14]. This model requires a description of opti- 
cal modes present in photonic crystals and dielectric micro- 
cavities. For example, in order to successfully describe the 
entanglement between photons and their quantum sources, it 
is imperative to achieve a resolution high enough to describe 
whispering gallery modes [15] available in dielectric micro- 
cavities such as micro-disks. This is especially important in 
applications that contrast classical computers, which depend 
on bits to store and process information, to quantum comput- 
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ers that depend on quantum bits (qubits). 

As the name suggests, quantum binary digits are a quantum 
representation of the on and off state as interpreted in machines 
like ENIAC 1 and those in existence today. It is possible to vi- 
sualize the relationship between bits and qubits by means of 
what is known as the Bloch sphere in terms of the on |1) and 
off states |0) of a bit and the state of a qubit One major 
advantage of quantum computers is that they do not alter the 
Church-Turing thesis [16] since the do not allow the compu- 
tation of functions which are not theoretically computable. So 
far it has been shown through the discoveries of Shor and oth- 
ers that it is possible to develop quantum algorithms for im- 
portant problems like prime factorization [17], protocols for 
quantum error correction (QEC), and fault-tolerant QC [18]. 
Other algorithms in QC, that if physically implemented could 
be of immediate use, include Grover's Algorithm for database 
searches [19] and the quantum Fourier Transform [20]. Also 
central to the discussion on QC and QEC is the decoherence 
rate of qubits. It is imperative to QC to find an implementation 
where qubits are well isolated from their environment [16]. 
Among suggested implementations for QC are Raman cou- 
pled low-energy states of trapped ions [3] and nuclear spins in 
chemical solution [4, 5]. Qubits based on these implemen- 
tations could provide the first examples of QC up to the 5 
through 10 qubits level. However, these implementations may 
not be scalable to more than 100 qubits [1]. 

Proposed implementations and promising schemes that 
could be scalable to more than 100 coupled qubits may be 
based on electron spins coupled by means of an optical mode 
of a photonic crystal or dielectric micro-cavity. One such 
scheme couples the electronic spin states of a Quantum Dot 
(QD) to the optical modes of a micro-disk [1]. Another 
promising scheme couples electronic QD states embedded 
inside nanocavities to the modes of a photonic crystal host 
[21, 22, 23]. An additional scheme recently realized experi- 
mentally has shown that Nitrogen- Vacancy centers in diamond 
can also be embedded inside a photonic crystal, which could 
enable fully scalable room-temperature quantum computing as 
a good alternative to using quantum dots. [24, 25, 26, 27]. 

The discussion for generating a theory on the interaction 
between such quantum sources and photonic states may be 
modeled after atomic systems [28]. For the case of a quan- 
tum source modeled after an effective two-level QD, the se- 
lection rules presented in figure (1) yield an interaction via the 
dipole approximation by implementing the state to state tran- 
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Fig. 1 : Selection rules in a two-level quantum source coupling 
to single photon states 



sition dipole moments of the QD [2] . To this end it is imper- 
ative to model spontaneous emission of a photon coupled to 
the electronic state of a quantum dot beyond the Weisskopf- 
Wigner approximation; which fixes the value of accessible 
modes from a possibly infinite set of frequencies to a single 
mode of a cavity or the transition frequency of the QD state as 
represented by uj 1 [28]. 

Our proposed model describes the photon by means of 
a Dirac-like equation for the photon. Experimentally, the 
quantum state of a photon may be reconstructed using op- 
tical homodyne tomography techniques by measuring quan- 
tum noise statistics of field amplitudes at different optical 
phases [29, 30]. In this work, the rigorous description of 
the interaction between a quantum source and the generated 
Maxwell Field is initially made within the formalism of rel- 
ativistic quantum field theory (QFT). In this description we 
begin within the canonical quantization procedure presented 
by the Gupta and Beuler method and the resulting interaction 
between these fields [31]. This procedure requires the defini- 
tion of a Lagrangian and gauge for the interacting fields 2 . To 
follow this procedure, a connection between the photon wave 
function (PWF) [32, 33, 34, 35] and the four vector potential 
for a Maxwell Field has to be drawn. This canonical quan- 
tization procedure leads to two important results, the com- 
plex Maxwell Field Tensor and the coupled electron-photon 
field equations in terms of a field equation for the PWF. We 
show the first result to be a Lagrangian for a free complex 
Maxwell Field written in terms of a self-dual tensor represent- 
ing the field tensor corresponding to the PWF that directly sat- 
isfies and yields equations of motion equivalent to the gener- 
alized Maxwell equations [36, 37, 38, 35]. In the second re- 
sult we show the coupling between a quantum source and the 
Maxwell Field it generates. We additionally extend these re- 
sults to show how these can be applied to model the emission 
of a single photon from a dielectric micro-cavity. Through- 
out the relativistic treatment of these fields we will maintain 
the Minkowski Metric to have the signature (+,—,—,—), and 
adopt the 4-notation consistent with x^ = (ct,x), r}^ v x v — 
x n (ct, x^ . 

Building on the work describing the PWF formalism, while 
working in Gaussian units in the presence of sources, one 



can define a self-dual tensor in terms of the electromagnetic 
field tensor [36] $^ = d^A v - d v and its dual = 

le^va^^ as W = 3^ - in terms of the vector P°" 
tential by 

0^ = - d v A») - l - e ^ (d,*A„ - d v A^) 



from which it is possible to define a general gauge-invariant 
Lagrangian for a free photon via 
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The same Lagrangian expressed in terms of the well known 
Faraday tensor [3 9] reads 
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with the corresponding matrix representation for the self-dual 
tensor given by 
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Using E and B to represent electric and magnetic fields re- 
spectively [40] , F T ee B ±iE [33, 32, 39]. This vector may 
be written in terms of the Riemann-Silberstein vector [34, 35] 
through the identity T± = =biF± = E ± iB (given that i 
represents the unit pseudo-scalar [41]). 

We defer to the traditional Quantum Electro-Dynamics 
(QED) interaction between an electron, a photon, and an addi- 
tional gauge field as mitigated by ^S^Tot^ where (Ijf. = 7 /i A /i ), 
via ^Tot = 4- + ^Ext and therefore express the interaction La- 



grangian for a gauge field A^ and spinor i/; 
to the the description in [42] by writing 
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It is worth noting that this Lagrangian leads to electromagnetic 
fields that satisfy the principle of superposition as required by 
experiment along with their conservation laws and definition 
of spin. This is evident because there are only expressions 
quadratic in the field and first order time derivatives present in 
the action. Explicitly in terms of the relativistic equations of 
motion and their Hermitian conjugates 
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In the relativistic regime we set V • A -\ — dt& = 0. 



In the non-relativistic limit working in the radiation gauge, 
these EOMs yield the well known Pauli-Schrodinger equation 
[43]. By quantizing the 3-vector potential while defining the 
new operators F+ V x A-i\d t ASc F _ VxA+i\d t A, 
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the EOM for the photon may be expressed as a Dirac-like 
equation 
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where a K k ' = —iSijk with representing the Levi-Civita 
permutation symbol. To define the quantized interaction term 
we expand the gauge field operator A in terms of creation and 
annihilation operators through the use of plane waves and re- 
tain a phase factor </> to account for its phase freedom [35] 
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Making the dipole approximation to the Pauli-Schrodinger 
equation [43], and changing to the interaction picture leads to 
the expression 
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Assuming that = c\k\ and making use of the identity [28] 

k 

— x x = —eric]: A , yields that the expression for the inter- 
action can be written in terms of the operator F+ 



V 



i(v k t-k-x) —i(<f>—% ) 



such that 



+ e* at e^t-k-x) e i{(j>-^)\ 

k,- k,- J 

by shifting the phase of the interaction « 

at xq, 

ihd t \(p) = - g J2 (^+ " • (pnmVnme^™ 1 ) \ip) (4) 

n,m 

The EOMs as derived from (4) (in terms of the photonic wave 
functions), for the case of a two level quantum source with an 
energy band-gap of AE = Hlj^, defined by the state- vector 3 
1^7) = Ca (t) \a0) + c b g (t) \blj*) interacting with it's own 
spontaneously emitted field, are given by 4 

ihc a (t)= ffil + ^) b e^ t -p ha 
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3 Where cr, 7 denote electronic and photonic states and a, b denote excited 
and ground states 
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Fig. 2 ! Analytically evaluated poles z n through use of Demoivre's theorem. 

Plots present the relative change of the argument 6 in arc-seconds 
A#" for zo & z\ . For Z2 the change in the argument is of the order 
A#" x 10 -2 . Radial length phasors z n + iu a for zo and z\ are 
of the order of 10 19 - and for zo these are of the order of 10 15 - 
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as plotted with respect to transition frequency uo a in 10 15 - and 
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with negative real components correspond to emission and those 
with positive real components correspond to revival in near field re- 
gions. 



Formally integrating these equations of motion while rep- 
resenting the polarization vectors as = (o + icj^j , 

eg _ = -±= (j) — icf^J , and the unit wave-vector as k = f, 
yield three poles z n which need to be considered when solv- 
ing for the probability density that ultimately gives rise to the 
spontaneously emitted photon, c b g A (£). Each of these poles 
correspond to emission and re-absorption respectively and al- 
low one to model Rabi oscillations associated with revival phe- 
nomena. These poles can be evaluated numerically or analyti- 
cally by application of Demoivre's theorem. 

In figure (2) we present solutions for poles z n in for transi- 
tion wavelengths and transition dipole moments in the ranges 
of 750 - 1300 nm and 20 to 100 Debye [44, 45, 46]. From 
these results, the poles of the contour integrals associated with 
the evaluation of c a (t) are presented as in figure (2) and their 
physical meaning interpreted from their location in the com- 
plex plane; where the contours are closed with further causal 
constraints. The expression for the wave-function, prior to its 
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Fig. 3i Near field revival phenomena represented in terms of energy ex- 
change AE (eV) between the QD state and the single photon for 
coarse and fine time steps. Envelope functions bounding the region 
of coherent oscillation behave as low order polynomials. The poly- 
nomial behavior of the envelope functions is contradictory to the 
expectation of exponential behavior with characteristic times of the 
order of the roots z n for c a . 



propagation through (3) is 



$*(+) 

7> — 5 o 



( 27 n=l 



( 27F ) n=l 

fto res^Zf we do not neglect terms of order O {r -2 }, usu- 
ally neglected in the far-field approximation [28]. Writing 
Pab,i -» pi along with Q n = uo a + zz n and * = x - x , 
and setting the speed of light within the interaction region to 
Co = r]~ 1 c, the components for an "emitting" pole with the 
contour closed over the lower half plane, are given by 
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as governed by the conditions that follow from the fact that the 
outgoing £_ and incoming £ + wave-fronts can not move faster 
than the speed of light. 
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We additionally go beyond the approximations which neglect 
revival [28], such that the components of an "absorbing" pole 
with the contour closed over the upper half plane, are given by 



# A-,)' 



(1±1) - 

~2~ Px 
(i±i) - 

o Px 



-£±,n,y — 



• ^ (lTl). 

Tlpx H Py 



±ipx 



2 

(1T1) 



Q 2 

tc 2 

n 2 



tci 



n 2 

^0 

n 2 



1 iQ n 
£ £ 2 c 



(p^ T Py) 
(px ± p y ) 



C4 



P^ 



^ =F Px) 
^ ± px) 





Px 


c+ 






2 




Px 


C- 


£c 


# - 


2 


iQ n 


Py' 


C+ 


£c 


# . 


2 


i£t n 


Py' 


C- 


£c 




2 



as governed again by the conditions that follow from the fact 
that the outgoing (_ and incoming £ + wave-fronts can not 
move faster than the speed of light 
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The coupling between (3), (5), and (6) was modeled compu- 
tationally through use of the following algorithm as imple- 
mented in a leap-frogging scheme [47] between real and imag- 
inary parts of the wave-functions, both within and beyond the 
interaction regions. 

• Initialize QD excited state 



• Use the leading coefficients and roots of the analytic form of c a 

& 



to determine the state of the photonic wave function ^ 
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at the next time step 
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• Repeat from step 2 
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Fig. 4i Finite difference models of the spontaneous emission of a single 
photon from a two-level system of quantum dot states. Color in- 
tensity in red (positive) and blue (negative) represent the amplitude 
of the photonic wave-function field strengths. 



Figure (3) shows the stable near-field exchange of energy be- 
tween the quantum dot and single photon states. These results 
are derived from two computational experiments set to have 
the same ratio of ^ . It is interesting to note that the envelope 
functions for both of these are very similar and both exhibit a 
seemingly linear decay of energy shortly after an initial revival 
of the quantum dot state. This seems to suggest that though 
it is worth while to retain a time resolution small enough to 
observe the initial revival event, this resolution can be made 
coarser almost immediately following the first revival. Figure 
(4) presents spatial XY and XZ plane projections of single 
photon state components F y ,+, F z + and their probability am- 
plitudes immediately following the establishment of coherent 
oscillations within the polynomial envelope functions shown 
in Figure (3). It is evident that during this time, though evanes- 
cent modes appear to escape the quantum dot region, calcula- 
tion of the localization of the single photon wave-function find 
the single photon state to be localized to the quantum dot re- 
gion in figures (4a) & (4b). The dipolar structure in the pro- 
jections presented by figures (4c) & (4e) was expected due to 
the initial orientation solely along the z-axis of the transition 
dipole moment of the two associated quantum dot levels. 

We have demonstrated that it is feasible to study near 
field single photon emission within dielectric structures by 
means of the Riemann-Silberstein wave-function beyond the 



Weisskopf-Wigner approximation; both analytically and com- 
putationally. It was further demonstrated that the locality of 
a photonic state could be well described during spontaneous 
emission while energy is injected and exchanged between both 
single photon and quantum dot states. Test cases were directly 
compared for different values of Ax & At . From these it 
was determined that the polynomial envelope functions for co- 
herent oscillations were in agreement within the initial revival 
period of the quantum dot excited state. Furthermore, the the- 
oretical approximations made accurately yield analytic and in- 
tuitive insight to the periodicity of the initial decay and revival 
phenomena present in the near field limit. This work therefore 
makes it feasible to computationally design photonic states to 
be emitted and detected by solid-state quantum dots embedded 
within dielectric structures and to compare them to experimen- 
tal results by means of their corresponding density matrix and 
Wigner functions. 
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